Mechanism of temperature dependence of the magnetic anisotropy energy in ultrathin 

Cobalt and Nickel films 
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Temperature dependent FMR-measurements of Ni and Co films are analysed using a microscopic 
theory for ultrathin metallic systems. The mechanism governing the temperature dependence of the 
magnetic anisotropy energy is identified and discussed. It is reduced with increasing temperature. 
This behavior is found to be solely caused by magnon excitations. 



I. INTRODUCTION 



Research on ultrathin films has been growing consid- 
erably over the last two decades due to their techni- 
cal importance and the increasing ability to grow high- 
quality film samples. When dealing with systems of re- 
duced dimensionality it is important to take into account 
the influence of magnetic anisotropics. Formally being a 
prerequisite for finite temperature magnetism 1 - they also 
play a crucial role in reorientation transitions as functions 
of film thickness, temperature, or an external magnetic 
field 2 .. 



In this work we present a microscopic quantum me- 
chanical model to describe the temperature dependence 
of the magnetic anisotropies of ultrathin Ni and Co 
films. As the dominant anisotropy terms in these tran- 
sition metals we take into account a second order lattice 
anisotropy and dipolar coupling. We are able to fit ex- 
perimental data obtained by the ferromagnetic resonance 
technique (FMR) 3 i4 with our model, which is based on 
a Heisenberg exchange between localized magnetic mo- 
ments. The main result will be that the T-dependence 
of the anisotropy is solely due to magnon excitations 
rather than to other mechanisms as, e.g., thermal ex- 
pansion or phononic interactions. The isolation of the 
magnon effect is of course not possible when describing 
the films with the classical (T=0-) Landau-Lifshitz equa- 
tions in which the (effective) anisotropy parameters have 
to be fitted at each given temperature^. Temperature is 
taken into account explicitly by the so-called stochastic 
Landau-Lifshitz-Gilbert (LLG) equation^. However we 
will introduce in the following an alternative approach 
based on a quantum mechanical description of the film 
system and an explicite consideration of magnon exci- 
tations using a Heisenberg model containing atomistic 
anisotropy terms. 



II. THEORETICAL DESCRIPTION 

The Hamiltonian of our microscopic model reads 
H = — JijSiSj — 2^ SJMsBoSj — 

ij i 
i ij \ *. 
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The first term describes Heisenberg coupling Jy between 
magnetic spin moments S; at the sites of a monolayer. 
It represents the largest energy scale in the problem and 
is responsible for the magnetism in the system. A film 
thickness beyond monolayer is effectively absorbed into 
the nearest neighbor exchange parameter J to which we 
restrict ourselves. J is chosen such that the monolayer 
magnetic moment equals that of the multilayer film at 
room temperature (T^ 7 = 410 K, T%° 2 = 400 K). The 
second term contains an external magnetic field Bo in ar- 
bitrary direction with the Lande factor gj and the Bohr 
magneton /is- The third and fourth term constitute lat- 
tice anisotropy and dipolar interaction, respectively, the 
latter leading to shape anisotropy. and g§ are micro- 
scopic anisotropy parameters, Si Z > is the z'-component of 
Si perpendicular to the film plane, and is the vector 
between lattice sites i and j . The shape anisotropy favors 
in-plane orientation and the lattice anisotropy can favor 
in-plane (K2 < 0) or out-of-plane (K2 > 0) orientation 
of the magnetization. 

The main idea of the method that we used in order 
to solve Q is discussed in Rcf.-. The exchange terms 
are decoupled using the standard Tyablikov (RPA) ap- 
proximation. The crucial point is to find a reasonable 
decoupling of the lattice anisotropy terms in the equa- 
tion of motion for the spin Green function ((S^~;S~)) 
given an arbitrarily oriented external magnetic field. This 
problem is solved by performing a coordinate transfor- 
mation (x',y',z r ) — ► (x,y,z). More precisely one self- 
consistcntly rotates the initial coordinate system defined 
by a z'-axis parallel to the film normal such that the z- 
axis of the new reference frame has the direction of the 
magnetization. Then an Anderson-Callen decoupling to 
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FIG. 1: Resonance frequency as a function of the external 
magnetic field applied parallel (6b q = 90°) and perpendicular 
(6b = 0°) to the film plane (monolayer). The resonance 
fields can be detected by an FMR experiment. Parameters: 
T = 0, K-2 = tO^skC, 50 = 0. 

the i^-terms is applied^. The results for the magneti- 
zation components obtained by this approach compare 
very well with corresponding QMC calculations 2 . We 
have improved the theoretical treatment in the mean- 
time by taking into account additional Green functions, 
namely all combinations of ((S' l + ' _ ; S^~'~)), as was pro- 
posed in Ref. ' . This improved theory ensures correct 
softening properties of the uniform spin wave mode and 
also agrees nicely with the QMC results. Furthermore 
the dipole term in (JJJ is treated in the RPA approxi- 
mation and only the uniform (q — > 0) contribution is 
considered as the non-uniform terms are negligible com- 
pared to contributions from the much larger Heisenberg 
exchange. 

A more detailed and general account of our method 
(e.g. the explicit extension to the multilayer case) will 
be presented elsewhere. We summarize only the essential 
output here: solving for the spin Green functions yields 
weights Xa(q) and excitation energies E a (q) = hui a (q) 
which in turn give the average magnon occupation num- 
ber 

q a 

The two terms of the sum over a describe the single- 
magnon excitations of the system for a given wave vec- 
tor q, namely magnon creation and magnon annihilation. 
The magnetization (in the rotated frame) can then be 
computed from 

(q v _ (1 + ^) 2S+1 (S - jg) + V 2S+1 (S + l+ jg) 

{ Z> (1 + ^)^+1-^25+1 • W 

Our theory therefore allows for a self-consistent determi- 
nation of the magnetization, i.e. temperature-dependent 
calculations, fn addition the self-consistent determina- 
tion of the rotation angle is achieved by requiring that 
S z be a constant of the motion, ^f- — 0, within the de- 
coupling approximations we use in our theory. We point 
out that by properly rescaling the parameters our equa- 
tions at T = can actually be shown to reduce to the 
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FIG. 2: Resonance field of Nickel at T = 55K and at T = 
294K as a function of the orientation of the external magnetic 
field. Circles: experimental data from Ni7/Cu(001). Lines 
connect theoretical fit points. S=l, Ki = 3.0^iBkG, go = 
4.5/xskG, J = 30meV. 



Landau-Lifshitz equations. A discussion of the impor- 
tant differences that appear for finite temperatures will 
be given elsewhere. 

In our theory the anisotropies K 2 and g influence the 
system solely via the effective anisotropy given by the 
temperature-dependent term 

K 2 (T) = (S Z )(T) (2K 2 C(T) - Dg ) , (4) 
C(T) = 1 _(^ + 1)- 2 (^)(T))^ (5) 

Here the T-independent quantity D is some number de- 
pending on the lattice geometry. Note that our effec- 
tive anisotropy K 2 is identical to the quantity M e ff 
commonly used within a Landau-Lifshitz description of 
FMR experiments 3 . Furthermore we exploit (S Z )(T) = 
S{S + 1) - (S Z )(T)(1 + 2<p{T)). The temperature de- 
pendence of K 2 thus goes beyond a mere proportionality 
to (S Z )(T) due to the occurence of the higher order In- 
dependent correlation function (S Z )(T). 

The experimental data have been obtained using FMR 
measurements 3 ^. This technique probes the uniform 
spin wave mode w(q = 0) of a magnetic sample. An ex- 
ternal field is tuned for a given probe frequency vq = 
o>(q=0)/27T until resonance occurs at B ics (9b ), with 
9b being the angle between the magnetic field and the 
normal to the filmplane. This is illustrated in Fig. ^ 

Using the temperature dependent effective anisotropy 
we can now fit the experimental data. We can 
thus check if the temperature dependence of the effec- 
tive anisotropy is due to spin wave excitations which are 
considered explicitly in our model or due to other 
(non-magnonic) effects which would manifest themselves 
in a temperature dependence of the parameters K 2 and 
9o- 
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- T = 55 K theoiy 
• T = 55 K exp. 

- - T = 294K theory 
o T = 294K exp. 



FIG. 3: Same as in Fig. 
Parameters: S'=2.5, K2 
J = 4.1meV. 
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2 for Cobalt (Cu 9 /Co 2 /Cu(001)). 
= -20.25^ s kG, go = 1.95/i B kG, 



III. FIT TO EXPERIMENTS WITH NI AND CO 
FILMS 

In the following we take S = 1 for Ni and S = 2.5 for 
Co due to different magnetic moments of the two metals^. 
The Lande factor is taken as gj = 2.1 for both Ni and 
Co. The FMR microwave frequency was set to 9 GHz. 

Fig. and 02 show the comparison between the 
B res (0B o )-cuTves from theory and experimental data 
for a (subscript denotes the number of monolayers) 
Cu 9 /Co 2 /Cu(001) and a Ni 7 /Cu(001) film system, re- 
spectively, at two different temperatures. In both cases 
the effective anisotropy favors the magnetization to lie 
within the film plane. There is quite good agreement at 
both temperatures over the whole range of angles 9b for 
both films. At a given angle the resonance field increases 
with temperature. 

It is important to note that the choice of the micro- 
scopic parameters K2 and go at a given temperature can- 
not be unambiguous as one easily sees from 10} . However 
we took additionally into account the temperature depen- 
dence of the resonance field for a fixed angle 9b = 90° 
as it is shown in Fig. 21 Due to the temperature depen- 
dent term which goes with K 2 in (gj, namely (S^)(T), 



the ambiguity is removed. Indeed it is still possible to 
accurately fit the experimental results with one set of (T- 
independent) parameters (if 2, go) for Ni and Co, respec- 
tively, over the whole temperature range. Furthermore in 
both cases the values of goS lie slightly above the result 
of an explicit evaluation of this quantity assuming point- 
like dipoles on the lattice sites for the given geometry^ 
(goS = 3.81/xskG). The conclusion we can draw is that 
the temperature dependence of the magnetic anisotropy 
energy is solely due to spin wave excitations which man- 
ifest themselves in the T-dependence of 0} rather than 
due to thermal expansion or phononic interactions. In 
other words, there is no additional T-dependence of the 
paramete rs K 2l go to be considered in order to des cribe 
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FIG. 4: Full temperature dependence of the resonance field 
in the easy direction (6b = 90°). Circles are experimental 
data, lines are from theory with the same parameters as in 
Fig. 2 and Fig. 3. 



the non-magnonic effects. 

In conclusion we presented a quantum theory for thin 
metallic films based on a local-moment model with lattice 
and shape anisotropy. By comparison with FMR experi- 
ments we found that the temperature dependence of the 
magnetic anisotropy energy of thin Ni and Co films is ex- 
clusively due to magnon excitations rather than caused 
by other structural or phononic effects. 

Financial support by the SFB 290 for part of this work 
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